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Abstract 



We have performed a large observing campaign of intermediate-redshift 
disk galaxies including multi-object spectroscopy with the FORS instru- 
ments of the Very Large Telescope and imaging with the Advanced Cam- 
era for Surveys onboard the Hubble Space Telescope. Our data set com- 
prises 113 late-type galaxies in the redshift range 0.1 < z < 1.0 and 
thereby probes galaxy evolution over more than half the age of the uni- 
verse. Spatially resolved rotation curves have been extracted and fit- 
ted with synthetic velocity fields that account for geometric distortions 
and blurring effects. With these models, the intrinsic maximum rota- 
tion velocity Knax was derived for 73 spirals within the field-of-view of 
the ACS images. Combined with the structural parameters from two- 
dimensional surface brightness profile fitting, the scaling relations (e.g., 
the Tully-Fisher Relation) at intermediate redshift were constructed. 
The evolution of these relations offers powerful tests of the predictions 
of simulations within the Cold Dark Matter hierarchical scenario. 

By comparing our sample to the Tully-Fisher Relation of local spi- 
ral galaxies, we find evidence for a differential luminosity evolution: the 
massive distant galaxies are of comparable luminosity as their present- 
day counterparts, while the distant low-mass spirals are brighter than 
locally by up to >2 m in rest-frame B. Numerous tests applied to the 
data confirm that this trend is unlikely to arise from an observational 
bias or systematic errors. Discrepancies between several previous stud- 
ies could be explained as a combination of selection effects and small 
number statistics on the basis of such a mass-dependent luminosity evo- 
lution. On the other hand, this evolution would be at variance with the 
predictions from numerical simulations. For a given V ma x, the disks of 
the distant galaxies are slightly smaller than those of their local coun- 
terparts, as expected for a hierarchical structure growth. Hence, the 
discrepancy between the observations and theoretical predictions is lim- 
ited to the properties of the stellar populations. A possible explanation 
could be the suppression of star formation in low-mass disks which is 
not yet properly implemented in models of galaxy evolution. 
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1 Introduction 



Within the last few years, our knowlegde of the basic parameters which deter- 
mine the past, present and future of the universe has improved significantly. 
Thanks to the combined results from studies of the Cosmic Microwave Back- 
grund, the Large Scale Structure, Big Bang Nucleosynthesis and distant su- 
pernovae, we now have strong evidence for a flat metric of spacetime (Spergel 
et al. 2003 and references therein). According to the observations, 73% of the 
mean density of the universe originate from Dark Energy, 23% are contributed 
by Cold Dark Matter and only 4% by "ordinary" baryonic matter. In such 
a cosmology, structure growth proceeds hierarchically, with small structures 
forming first in the early cosmic stages, followed by the successive build-up of 
larger structures via merger and accretion events. 

Although the constituents of the Dark Energy and Dark Matter remain 
unknown, the ACDM or "concordance" cosmology has been a very successful 
tool for the interpretation of structures on Mpc scales and beyond (e.g., Pea- 
cock 2003). On scales of individual galaxies, however, several discrepancies 
between observational results and theoretical predictions have been found, a 
prominent of which is the "angular momentum problem". This term depicts 
the loss of angular momentum of the baryons to the surrounding DM halo, re- 
sulting in galactic disks within numerical simulations which are smaller than 
observed (e.g., Navarro & White 1994), however more recent studies made 
progress in this respect (e.g., Governato et al. 2004). Aiming at a quanti- 
tative test of the hierarchical scenario at the scale of individual galaxies, we 
performed an observational study which covers a significant fraction of the 
Hubble time. 

For this purpose, we utilised scaling relations like the Tully-Fisher rela- 
tion (TFR, Tully & Fisher 1977) between the luminosity L and the maximum 
rotation velocity V maK of spiral galaxies. Basically, this correlation can be 
understood as a combination of the virial theorem and the rotational stabili- 
sation of late-type galaxies. By comparing local and distant spirals of a given 
Vmax, the luminosity evolution within the look-back time can be determined. 
Since the maximum rotation velocity is a measure for the total (virial) mass of 
a disk galaxy (V^ ax oc M v - lr , e.g. van den Bosch 2002), the TF analysis relates 
the evolution of stellar population properties to the depth of the gravitational 
potential well. 

Numerical simulations within CDM-dominated cosmologies have been suc- 
cessfully used to reproduce the observed slope of the local TFR, whereas the 
zero points were offset due to dark halos with too high concentrations (e.g., 
Steinmetz & Navarro 1999). The TFR slope was predicted to remain con- 
stant with cosmic look-back time in such A-body simulations; nevertheless 
the modelling of realistic stellar populations at sufficient resolution remains 
a challenge — typically, the masses of individual particles are of the order of 
10 5 Mq . . . 10 6 M Q . Other theoretical approaches focussed more on the chemo- 
spectrophotometric aspects of disk galaxy evolution. For example, Boissier 
& Prantzos (2001) calibrated their models to reproduce the observed colors 
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of local spirals. Compared to these, the authors predicted higher luminosi- 
ties for massive disks and lower luminosites for low-mass disks at redshifts 
z > 0.4. A similiar evolution was found by Ferreras & Silk (2001). By mod- 
elling the mass-dependent chemical enrichment history of disk galaxies with 
the local TFR as a constraint, the authors found a TFR slope that increases 
with look-back time (i.e., for a parameterisation L oc V^axi a increases with 
redshift) . 

In the last decade, many observational studies of the local TFR have pro- 
duced very large samples with iV Q bj « 1000 (e.g. Haynes et al. 1999), not 
only to derive the slope and scatter with high accuracy, but also to map the 
peculiar velocity field out to cz w 15000 km s _1 (e.g. Mathewson & Ford 
1996). Other groups used spirals, partly with cepheid-calibrated distances, 
to measure the Hubble constant. For example, Sakai et al. (2000) derived a 
value of H = (71 ± 4) kms -1 Mpc -1 with this method. 

At higher redshifts, robust measurements of rotation velocities become 
increasingly difficult, which is mainly for two reasons. Firstly, because the 
objects are very faint. Given a redshift of z — 0.5, the surface brightness 
at galactocentric radii of ~ 3rd — where the regime of constant rotation 
velocity is reached — is typically /is ~ 27 mag arcsec" 2 . Spatially resolved 
spectroscopy at this level has become feasible just with the generation of 10m- 
class telescopes. The second difficulty coming into play arises from the small 
apparent sizes of the galaxies, this issue will be described in detail in Sect. 21 

A number of samples with 10-20 objects in the regime 0.25 < (z) < 0.5 
have been observed to estimate a possible evolution in luminosity by com- 
parison to the local TFR. The results of these studies were quite discrepant: 
e.g. Vogt et al. (1996, 1997) found only a modest increase in luminosity of 
AM B w -0.5 m , whereas Simard & Pritchet (1998) and Rix et al. (1997) 
derived a much stronger brightening with AMb ~ — 2.0 m . A more recent 
study of 19 field spirals by Milvang-Jensen et al. (2003) yielded a value of 
AMb ~ —0.5™ and showed evidence for an increase of this brightening with 
redshift. 

It seems likely that some of these results are affected by the selection cri- 
teria. For example, Rix et al. selected blue colors with (B — i?) b s < 1.2 m , 
Simard & Pritchet strong [Oil] emission with equivalent widths >20 A, while 
Vogt et al. partly chose large disks with > 3 kpc. The two former criteria 
prefer late-type spirals, whereas the latter criterion leads to the overrepresen- 
tation of large, early-type spirals. Additionally, due to the small samples, all 
these studies had to assume that the local TFR slope holds valid at interme- 
diate redshifts — we will adress this topic again in Sect. HO 

Based on a larger data set from the DEEP Groth Strip Survey (Koo 2001) 
with N w 100 spirals in the range 0.2 < z < 1.3, Vogt (2001) found 
a constant TFR slope and only a very small rest-frame B-band brightening 
of < 0.2 m . On the other hand, in a more recent study based on the same 
survey, an evolution of the luminosity-metallicity relation both in slope and 
zero point was observed (Kobulnicky et al. 2003). The authors argued that 
low-luminosity galaxies probably have undergone a decrease in luminosity 
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combined with an increase in metallicity during the last ~ 8 Gyr. 

Throughout this article, we will assume the concordance cosmology with 
n m = 0.3, A = 0.7 and H = 70kms~ :L Mpc" 1 . 

2 Sample Selection & Observations 

The sample described here has been selected within the FORS Deep Field 
(FDF, see Heidt et al. 2003), an UBgRIJK photometric survey covering 
a sky area of ~ 6 x 6 arcmin 2 near the southern Galactic pole. The imag- 
ing was performed with the Very Large Telescope (optical bands) and the 
New Technology Telescope (Near Infrared bands) . Based on a catalogue with 
spectral types and photometric redshift estimates (Bender et al. 2001), we 
chose objects for follow-up spectroscopy which satisfied the following crite- 
ria: 1) late-type Spectral Energy Distribution, i.e., galaxies with emission 
lines, 2) total apparent i?-band magnitude R < 23 m , 3) photometric redshift 
z P hot < 1-2 to ensure that at least the [O Il]3727 doublet falls within the wave- 
length range of the spectra, 4) disk inclination angle i > 40° and 5) deviation 
between slit direction and apparent major axis of 6 < 15°. The two latter 
constraints were chosen to limit the geometric distortions of the observed ro- 
tation curves. For some objects, however, these limits had to be exceeded 
during the construction of the spectroscopic setups. 

After a pilot observation in 1999, the spectroscopy was performed in 2000 
and 2001 using the FORS1 & 2 instruments of the VLT in multi-object spec- 
troscopy mode with a total integration time of 2.5 hrs per setup. Using the 
grism 600R, a spectral resolution of R s=s 1200 was achieved with a spectral 
scale of 1.07 A/pix and a spatial scale of 0.2 arcsec/pix. The seeing ranged be- 
tween 0.4 and l.Oarcsec with a median of 0.74arcsec. In total, 129 late-type 
galaxies were observed. 

For an accurate derivation of the galaxies' structural parameters, like disk 
inclination, scale length etc., we also took Hubble Space Telescope images of 
the FDF with the Advanced Camera for Surveys using the F814W filter. To 
cover the complete FDF area, a 2 x 2 mosaic was observed. 

3 Spectrophotometric Analysis 

The spectra of 113 galaxies were reliable for redshift determination. Out of 
these, 73 objects eventually yielded maximum rotation velocities (see next 
section) and were covered by the HST/ACS imaging; these objects will be 
referred to as the FDFTF sample in the following. They span the redshift 
range 0.09 < z < 0.97 with a median of (z) — 0.45 corresponding to look-back 
times 1.2 Gyr < t\ < 7.6 Gyr with (ti) = 4.7 Gyr. This data set covers all 
spectrophotometric types from very early-type spirals (Sa or T — 1) to very 
late-type galaxies (Sdm/Im or 8 < T < 10). 

An analysis of the galaxies' surface brightness profile profiles was con- 
ducted with the GALFIT package (Peng et al. 2002). To fit the disk com- 
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ponent, an exponential profile was used, while a potential bulge was approxi- 
mated with a Sersic profile. In the case of 13 FDFTF galaxies, the fit residual 
images and large fit errors indicated an irregular component that could not 
be approximated properly with a Sersic law. The bulge-to-total ratios of 
these galaxies were assumed to be undefined. The B/T ratios of the other 60 
FDFTF galaxies (0 < B/T < 0.53 with (B/T) = 0.04) confirm that the vast 
majority of these galaxies are disk-dominated. 

Total apparent magnitudes were determined using the mag_auto algorithm 
of the Source Extractor package (Bertin & Arnouts 1996). For the computa- 
tion of absolute -B-band magnitudes Mb, we used the filter which, depending 
on the redshift of a given object, best matched the rest-frame -B-band. For 
galaxies at z < 0.25, 0.25 < z < 0.55, 0.55 < z < 0.85 and z > 0.85, we thus 
utilised the B, g, R and I magnitudes, respectively. Thanks to this strategy, 
the fc-correction uncertainties Ok — usually a substantial source of error to 
the luminosities of distant galaxies — are smaller than 0.1™ for all types and 
redshifts in our sample. For the correction of intrinsic dust absorption, we 
followed the approach of Tully & Fouque (1985) assuming a face-on (i = 0°) 
extinction of As = 0.27" 1 . The absolute magnitudes of the FDFTF galaxies 
computed this way span the range — 18. m > Mb > — 22. 7 m . 

The spectra of 12 objects in our sample cover a wavelength range that 
simultaneously shows emission in [Oll]3727, H/3, [OlIl]4959 and [Olll]5007 
at sufficient signal-to-noise to determine the equivalent widths. These lines 
can be used to estimate the gas-phase metallicity. We adopted the analytical 
expressions given by McGaugh (1991) to compute the abundances O/H from 
the i?23 and O32 parameters. Since all the galaxies have Mb < —18™, we 
assumed that they fall on the metal- rich branch of the -R23-O/H relation. 
The galaxies have abundances 8.37 < log(0/H) < 8.94. We will use these 
estimates to investigate the luminosity-metallicity relation in Sect.|3 



4 Derivation of V max 

We extracted spatially resolved rotation curves from the two-dimensional 
spectra by fitting Gaussians to the emission lines stepwise along the spatial 
axis. Line fits at any projected radius which, compared to the instrumental 
broadening (FWHM ins w 4.5 A), had very small (FWHM fit < 2 A) or very 
large (FWHM fit > 12 A) line widths were assumed to be noise and therefore 
neglected. 

The analysis of spatially resolved rotation curves from optical spectroscopy 
of local spiral galaxies is relatively straightforward. But in the case of distant 
galaxies with very small apparent sizes, the effect of the slit width on the 
observed rotation velocities V I0 t(r) must be considered. At redshift z = 0.5, a 
scale length of 3 kpc — typical for an L* spiral — corresponds to ~ 0.5 arcsec 
only, which is half the slit width used in our observations. Any value of V I0 t(r) 
is therefore an integration perpendicular to the spatial axis (slit direction), 
a phenomenon which is the optical equivalent to "beam smearing" in radio 
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Figure 1: Example of a simulated rotation velocity field for an object from our 
data set with a disk inclination i = 64° and a misalignment angle between 
apparent major axis and slit direction of S = +13°. The dashed ellipses 
denote the isophotes of the disk, with the outermost corresponding to an I- 
band surface brightness of [ii « 25magarcsec~ 2 . The curved dotted and solid 
lines correspond to line-of-sight rotation velocities ranging from — 120km/s 
to +120km/s. The two solid horizontal lines visualise the position of the slit 
used for spectroscopy. 



observations. The seeing has an additional blurring effect on the observed 
rotation curves. If not taken into account, these two phenomena would lead 
to an underestimation of the intrinsic rotation velocities and, in particular, 
the intrinsic V max . 

We overcame this problem by generating synthetic rotation curves. For 
the intrinsic rotational law, we used a simple shape with a linear rise of V TO t (r) 
at small radii, turning over into a region of constant rotation velocity where 
the Dark Matter halo dominates the mass distribution. Alternatively, we also 
tested the so-called "Universal Rotation Curve" shape (Persic et al. 1996), a 
parameterisation which introduces a velocity gradient in the outer regions of 
the disk which is positive for sub-L* objects and negative for objects much 
more luminous than L*. However, the results given here are not sensitive to 
the form of the intrinsic rotational law — see Bohm et al. (2004) for a detailed 
discussion of this topic — and we therefore only use V max values determined 
with an intrinsic "rise-turnover-flat" shape here. 

Given the observed inclination, position angle and scale length of an object, 
the intrinsic rotation velocity field was constructed, Fig. shows an example. 
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Figure 2: Examples of rotation curves from our data set. The solid lines 
are the synthetic rotation curves fitted to the observed rotation velocity as a 
function of radius (solid symbols) used to derive the intrinsic maximum rota- 
tion velocity. For each object the spectrophotometric type, total apparent R 
magnitude, redshift and Vmax are given. The two upper curves were classified 
as high quality data, the two lower ones as low quality data due to the large 
measurement errors (lower left) or an asymmetric shape (lower right). 



In the next step of the simulation, the velocity field was weighted with the 
surface brightness profile. The effect of this was that, just like for the observed 
data, brighter regions contributed stronger to the rotation velocities in direc- 
tion of dispersion than fainter regions (the "beam smearing" effect) . Following 
the weighting, the velocity field was convolved with the Point Spread Func- 
tion to simulate the blurring due to seeing. Finally, a "stripe" was extracted 
from the velocity field, with a position and width that corresponded to the 
slit used during the observations, and integrated perpendicular to the spatial 
axis. The results of the whole procedure was a synthetic rotation curve which 
introduced the same geometric and blurring effects as the corresponding ob- 
served rotation curve. By fitting the simulated rotation curve to the observed 
rotation curve, we derived the intrinsic value of Vmax. Four examples of ob- 
served rotation curves along with the best-fitting synthetic rotation curves are 
shown in Fig. 

36 galaxies had to be rejected from the further analysis because the S/N 
was too low to probe the regime of constant rotation velocity at large radii, 
or because the rotation curves were perturbated. Four objects were reliable 
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Figure 3: FORS Deep Field sample of spirals in the range 0.1 < z < 1.0 in 
comparison to the local Tully-Fisher relation by Pierce & Tully (1992); the 
solid line denotes the fit to the local data, the dashed lines give the 3 a limits. 
The distant sample is subdivided according to rotation curve quality: high 
quality curves (solid symbols) extend well into the region of constant rotation 
velocity at large radii and therefore give robust values of V max , Error bars are 
shown for the high quality data only. 



for the Vmax determination, but their positions were located outside the field 
of-view of the HST/ACS mosaic imaging. 34 objects had curves with a high 
degree of symmetry and clearly reached into the "flat" regime, we consider 
these as high quality data. 39 curves had a relatively small spatial extent or 
mild asymmetries, these will be referred to as low quality data in the following. 
In total, our kinematic data set thus comprises 73 late-type galaxies at a mean 
look-back time of ^5Gyr. The objects span the range 25 km/s < V max < 
450 km/s with a median of 129 km/s (high quality data only: 62 km/s < 
Vmax < 410 km/s and (V max ) = 154 km/s). 



5 Scaling Relations at Intermediate Redshift 

In Fig. El the maximum rotation velocities and absolute magnitudes of the dis- 
tant FDFTF galaxies are compared to the local B-band Tully-Fisher relation 
by Pierce & Tully (1992): 

Mb = -7.48 log V max - 3.52 (1) 
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with a scatter of cr B = 0.41™. Note that, at variance with the original rela- 
tion given by these authors, we have calibrated the zero point to a face-on 
extinction of 0.27™ 1 to achieve consistency with the computation of the distant 
galaxies' absolute magnitudes. We emphasize that the further analysis is not 
sensitive to the choice of the local reference sample: e.g., for the large data 
set of Haynes et al. (1999, comprising 1097 objects), we find a very similar 
relation of 

M B = -7.85 log V max - 2.78, (2) 

using a bisector fit (two geometrically combined least-square fits with the 
dependent and indepedent variable interchanged). We will utilise the Pierce 
& Tully sample here for the sake of comparability to intermediate-redshift 
TF studies in the literature which mostly have used this sample as a local 
reference. 

On the average, the distant galaxies are overluminous with respect to their 
local counterparts, we find a median offset of (AMb) — —0.98™ for the total 
FDFTF sample and (AMb) = -0.81 m for the high quality data only. But 
we find also evidence for a differential evolution. Fig. indicates a relatively 
good agreement between the intermediate-redshift galaxies and the local TFR 
in the regime of fast rotators, i.e. high masses, while the distant low-mass 
galaxies systematically deviate from the relation of present-day spirals. For 
low quality data, this may partly be due to underestimated maximum rotation 
velocities, since the corresponding curves have a relatively small spatial extent 
and do not robustly probe the region of constant rotation velocity at large 
radii. In the case of high quality rotation curves, this is however unlikely, 
since these extent well into the "flat region". 

A 100 iteration bootstrap bisector fit (average of 100 bisector fits with 
randomly removed objects in each iteration) to the 34 FDFTF objects with 
high quality rotation curves yields 

M B = -(4.05 ±0.58) log ^ max - (11.8 ± 1.28), (3) 

i.e. the TFR slope we find at intermediate redshift is significantly shallower 
than in the local universe. Since the derivation of the galaxies' structural 
parameters and of the Knax values has been based entirely on HST/ACS 
imaging, Eq.0is a confirmation of the results presented in Bohm et al. (2004) 
which were limited to ground-based imaging. 

We show the individual offsets AMb of the FDF galaxies from the local 
TFR as a function of their maximum rotation velocity in Fig. Even when 
restricting the sample to the high quality rotation curves, we find significant 
overluminosties of up to more than 2™ 1 in the rest-frame B for low-mass spi- 
rals. L* galaxies, corresponding to log Vm ax ~ 2.3 according to Eq. scatter 
around a negligible evolution, while the fastest rotators are systematically 
under luminous. 

In Fig. El we show the sub-sample of 12 FDF galaxies for which we could 
determine the oxygen abundances O/H in comparison to the local luminosity- 
metallicity relation as given in Kobulnicky et al. (2003, the displayed scatter 
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Figure 4: Offsets of the distant FORS Deep Field galaxies from the local 
TFR by Pierce & Tully (1992). The distant sample is subdivided according 
to rotation curve quality: high quality curves (solid symbols) extend well into 
the region of constant rotation velocity at large radii and therefore give robust 
values of V max . The dashed horizontal line corresponds a zero luminosity 
evolution. While high-mass galaxies are in agreement with the local TFR 
or even slightly underluminous given their V max , the objects are increasingly 
overluminous towards small values of V max (error bars are shown for the high 
quality data only). 

is a rough estimate). In addition, a sub-sample of distant late-type galax- 
ies from the DEEP survey (ibid.) is shown, which has been restricted to 
the same redshift intervall (0.22 < z < 0.46) that is covered by the FDF 
galaxies. Both sub-samples thus represent a look-back time of ~ 4 Gyr. For 
the sake of comparability, the absolute magnitudes Mb of the FDF galax- 
ies given in this figure are not corrected for intrinsic absorption, as is the 
case for the Kobulnicky et al. data. Both distant samples indicate a "tilt" 
with respect to the local L-Z relation. At given log(0/H), high-metallicity 
galaxies at intermediate redshift agree relatively well with the local L-Z re- 
lation, whereas low-metallicity objects are overluminous. Alternatively, the 
distributions may be interpreted such that the distant low-luminosity galax- 
ies have smaller chemical yields than locally, while high-luminosity galaxies 
do not differ strongly in O/H between intermediate and low redshift. If the 
offsets of the FDF spirals we observe in the TF diagram are due to younger 
stellar populations than locally, it is probable that Fig. shows a combined 
evolution in luminosity and metallicity. This indeed has been the conclusion 
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Figure 5: The gas-phase metallicities of 12 spirals from the FDF data set 
(filled symbols) in comparison to the local luminosity-metallicity relation con- 
structed by Kobulnicky et al. (2003, solid line; dashed lines give the esti- 
mated la scatter). Also shown are distant galaxies presented in Kobulnicky 
et al. (open symbols) which cover a similar redshift range as the FDF galax- 
ies. Both distant samples show a "tilt" with respect to the local relation 
which likelywise indicates a combined evolution in luminosity and metallicity 
of low-luminosity galaxies. 

of Kobulnicky et al. after comparison of their data to single-zone models. 

The third scaling relation which we want to focus on here is presented in 
Fig. E3 where the FDF galaxies are shown with respect to the local velocity- 
size relation correlating Vmax and rd. To derive the latter, we used the sample 
of Haynes et al. (1999). A bisector fit this data set yields 

logr d =1.35 log V max - 2.41. (4) 

Since the scale lengths of the FDF galaxies were determined in the /-band 
(HST-filter F814W), the data intrinsically probe shorter wavelengths towards 
higher redshifts. For direct comparability to the local sample (/-band as well), 
we had to account for this rest-frame shift. Otherwise, the measured FDF 
disk sizes could be overestimated, in particular for the more distant galax- 
ies. Adopting the relations between scale lengths of local spirals at different 
wavelengths presented in de Jong (1996), we have transformed the observer's 
frame /-band scale lengths to the rest-frame /-band values. Note, however, 
that these factors are relatively small: for the FDF galaxies at z « 1 , the cor- 
rection corresponds to only ~ 10%. The characteristic disk sizes of the FDFTF 
sample cover the range 0.7kpc < rd < 10.1 kpc with a median (rd) = 2.7kpc. 
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Figure 6: Comparison between the intermediate-redshift FDF galaxies and 
the local velocity-size relation of the Haynes et al. sample (1999, solid line; 
dashed lines correspond to the 3 a scatter) . Error bars are shown for high 
quality data only. 

6 Discussion 

The Vmax-dependent TF offsets we observe at redshift (z) w 0.5 may be in- 
dicative for a significant decrease of the luminosity of low-mass galaxies — 
possibly combined with an increase in metallicity — over the past ~ 5 Gyr and 
a negligible evolution of high-mass galaxies. This evolution would be at vari- 
ance with theoretical predictions: e.g., Steinmetz & Navarro (1999) find mass- 
independent TF offsets towards higher redshifts with an A^-body Smoothed 
Particle Hydrodynamics code. Boissier & Prantzos (2001), who used a "back- 
wards approach" model calibrated to the observed chemo-spectrophotometric 
properties of local spirals, predict overluminosties of high-mass spirals and 
under luminosities of low-mass spirals towards larger look-back times. 

It has to be ruled out that our result might be induced by an observational 
bias or a systematic error. E.g., it is known that present-day spirals show a 
correlation between their TF residuals and broad-band colors (e.g. Kannap- 
pan et al. 2002), with blue galaxies preferentially populating the regime of 
over luminosities. We have therefore tested whether Fig. 31 may simply reflect 
an evolution of the color-residual relation with redshift, finding no evidence 
for such a trend (Bohm et al. 2004). Another issue that has to be adressed is 
the potential impact of sample incompleteness. Any magnitude-limited data 
set contains only a fraction of the objects that are located within the observed 
volume. Towards lower luminosities (or slower rotation velocities), this frac- 
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tion becomes smaller. Furthermore, the magnitude limit corresponds to higher 
luminosities at higher redshifts. An incompleteness bias could therefore result 
in a flattening of the distant TFR with increasing redshift. However, dividing 
our sample into objects with z < 0.449 (37 galaxies) and z > 0.449 (36 galax- 
ies), we find no evidence for such a redshift dependence, the respective slopes 
of the two redshifts bins are —3.49 and —3.77. For a more sophisticated test of 
sample incompleteness which is based on the work of Giovanelli et al. (1997), 
we refer to our results presented in Ziegler et al. (2002). 

In the following, we will address another three examples of tests we per- 
formed. These are related to the influence of the intrinsic rotation curve shape, 
the impact of the intrinsic absorption correction and the issue of galaxy-galaxy 
interactions. 

To derive the intrinsic maximum rotation velocity, we have assumed an in- 
trinsic rotation curve shape with a linear rise of the rotation velocity at small 
radii which turns over into a region of constant rotation velocity at a radius 
that depends on the rest-frame wavelength of the used emission line. This 
shape is observed for kinematically unperturbated, massive (~ L*) local spi- 
rals (e.g., Sofue & Rubin 2001). For galaxies of very high or very low masses, 
on the other hand, it is observed that even in the outer parts, most rota- 
tion curves have a velocity gradient. While the rotation velocity keeps rising 
beyond the "turnover" radius in very low-mass spirals, the velocity gradient 
in very high-mass spirals is negative. Persic et al. (1996) have used >1000 
curves of local spirals to derive a parameterisation that uses the luminosity 
of an object as an indicator for the rotation curve shape. To ensure that the 
observed TF offsets cannot be attributed to a false assumption on the intrinsic 
rotational law, we have alternatively used this so-called "Universal Rotation 
Curve" shape as input for the computation of our synthetic velocity fields. If 
we use the V maK values derived this way to recompute the offsets from the 
local TFR, the luminosity evolution we find is smaller by only ~0.15 m at 
V max ~ lOOkm/s. Since this is a very modest change of the offsets we found 
on the basis of the simple "rise-turnover-flat"-shape (which have a median of 
(AM B ) « -1.74 m at V max « lOOkm/s for the HQ data), we conclude that 
our results do not differ significantly between these two assumptions on the 
intrinsic rotation curve shape. 

Similarly, we have tested whether a different approach to correct for the 
intrinsic absorption would have an effect on our results. All values given here 
were derived following Tully & Fouque (1985), i.e., the amount of intrinsic 
absorption is assumed to depend only on the inclination of the disk. More 
recently, Tully et al. (1998) have found evidence that the dust reddening is - at 
least locally — stronger in high-mass spirals than in low-mass spirals. Using 
their results, we have recomputed the absolute magnitudes of the FDFTF 
galaxies. As a new local reference that is consistently corrected for intrinsic 
absorption following Tully et al., we adopted the sample of Verheijen (2001) 
which is slightly steeper (slope —8.1) than the Pierce & Tully (1992) sample. 
This is simply due to the fact that fast rotating, high-mass spirals are assumed 
to have a larger amount of intrinsic absorption than in the Tully & Fouque 
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approach, and vice versa in the low-mass regime. The offsets of the FDF high 
quality data from the Verheijen TFR are however very similar to the initial 
values ((AM B ) « -1.77" 1 vs. (AM B ) « -1.74 m at y max w lOOkm/s). With 
respect to the two conventions of intrinsic absorption correction discussed 
here, the TF offsets therefore are robust. 

A third aspect we want to focus on here concerns the interplay between 
the stellar population properties and the environment. From studies in the 
local universe, it is known that galaxies residing in close pairs can be subject 
to tidal interactions which can increase the star formation rates. In such 
cases, the fraction of high-mass stars would be enlarged and, in turn, the 
mass-to-light ratio would be decreased, resulting in overluminosities in the 
TF diagram. This triggering of star formation would be particularly efficient 
in low-mass galaxies (e.g., Lambas et al. 2003). Though we have selected 
our objects from a sky region that should be representative for low-density 
environments, it is not clear a priori whether the correlation between the TF 
offsets and V max can at least in part be attributed to tidally induced star 
formation. Based on >10 5 galaxies from the 2dF survey, Lambas et al. have 
found that the star formation rates can be significantly increased in objects 
that have close companions with a separation A14 ys < 250 km/s in systematic 
velocity and a projected distance of -D pro j < 100 kpc. Using all 267 available 
spectroscopic redshifts of FDF galaxies at z < 1 from our own study and 
a low-resolution survey presented in Noll et al. (2004), and adopting the 
Lambas et al. constraints cited above, we have found 12 FDFTF objects with 
confirmed neighbors. 

In Fig. we show the TF offsets of these galaxies in comparison to the 
rest of the sample. Though the small sub-sample of pair candidates does 
not allow robust statistics, the galaxies with close companions appear to be 
similarly distributed as the rest of the sample. Moreover, we find a hint that 
the rotation curve quality is reduced with respect to the probably isolated 
galaxies — only 3 out of 12 (25%) pair candidates were classified to have 
high quality rotation curves, whereas for the rest of the sample, this fraction 
is 31 out of 61 (51%). In particular, the pair candidates with HQ curves 
are not systematically biased towards large overluminosities. Since we have 
spectroscopic redshifts only for a fraction of the galaxies within the probed 
volume, it is possible that we missed some close pairs. However, the aim was 
to test whether those pair candidates which are identified systematically differ 
from the remaining FDFTF objects. We thus conclude that our analysis which 
is based on the high quality rotation curves is very unlikely to be affected by 
tidally induced star formation. 

To summarise all the tests performed, we find no evidence for any system- 
atic error or bias that may be the source of the observed shallow slope of the 
intermediate-redshift TFR. We now will show that our findings can be used 
to interpret the rather discrepant results of previous TF studies of distant 
galaxies introduced in Sect.Q] For this, Fig. shows our sample sub-divided 
according to the Spectral Energy Distribution into early-type spirals (Sa/Sb), 
intermediate-types (Sc) and very late-types (Sdm/Im). All three sub-samples 
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Figure 7: Offsets of the distant FORS Deep Field galaxies from the local 
TFR by Pierce & Tully (1992). Large symbols depict FDF objects which show 
spectroscopically confirmed neighbors within AV^ ys < 250 km/s in systematic 
velocity and a projected separation -D pro j < 100 kpc. See text for details. 



show a correlation between the TF offsets and V m&K , but cover different 
mass regimes: galaxies with late-type spectra have smaller average V max val- 
ues than early-type spirals. The respective classes have median values of 
(V r ma x)sdm/im = 91 km/s, (V r max )sc = 140 km/s and (V r max > Sa/sb = 240 km/s. 
For samples which are too small to robustly test a correlation between AMb 
and Vmax but can only be used to determine average TF offsets, this would 
results in a correlation between SED type and (AMb)- 

Studies comprising only 10-20 galaxies which were selected on blue colors 
or strong emission lines and therefore preferentially contained late-type spirals 
(e.g., Rix et al. 1997, Simard & Pritchet 1998), found evidence for large 
luminosity offsets at intermediate redshift with respect to the local TFR. 
According to the above, this can be attributed to a combination of a selection 
effect and small number statistics. Similary, a study by Vogt et al. (1996) 
which was selected on large disks and therefore mainly contained early-type 
spirals (which on average are more massive than late-types) , yielded a modest 
value for the mean TF offset. 

A more recent study on ^100 intermediate-redshift spirals by Vogt (2001) 
did not find evidence for a slope evolution of the TFR, at variance with 
our results. Moreover, this study yields only a modest average TF offset of 
— 0.2 m , compared to the — 0.8 m we find for the HQ data (note that these 
values correspond to the same cosmology and similar redshifts) . It is difficult 
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Figure 8: Offsets of the distant FORS Deep Field galaxies from the local 
TFR by Pierce & Tully (1992). The distant sample is subdivided according 
to the SED type. All three sub-samples indicate a AMB-V mavi correlation, 
but cover different mass regimes. See text for details. 



to speculate whether this could be attributed to different criteria of rotation 
curve quality or differences in the Vn iax derivation procedures. On the other 
hand, a sample of 64 galaxies drawn from the same survey (the DEEP project, 
see Koo 2001) finds, in comparison to local galaxies, a tilt of the intermediate- 
redshift luminosity-metallicity relation which is also indicated by our data 

(Fig. ED. 

Is is a complicated issue to determine the key processes which can give rise 
to the mass-dependent TF offsets we observe. Several effects likelywise act 
in combination when local and distant spirals are compared. E.g., the stellar 
populations of the intermediate-redshift galaxies are probably younger than 
those of their local counterparts, the gas mass fractions and chemical yields 
also evolve with time etc. Since even a relatively small fraction of young, 
high-mass stars can have a significant effect on the luminosity in the blue 
bands, a straightforward interpretation would be that the TF offsets of the 
FDF galaxies point towards a correlation between mass and age. 

As a first deeper insight into this matter, Ferreras et al. (2004) have used 
single-zone models of chemical enrichment on a sub-sample of the FDFTF 
galaxies at z > 0.5. These models were determined by only four free parame- 
ters: formation redshift, gas infall timescale, star formation efficiency and gas 
outflow fraction. Model star formation histories were generated which, com- 
bined with the latest Bruzual & Chariot models (2003), were used to compute 
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Figure 9: Offsets Alogrd of the distant FORS Deep Field galaxies (high 
quality rotation curves only) from the local velocity-size relation of the sample 
by Haynes et al. (1999) shown in Fig. plotted as a function of redshift. 
Objects with A logrj > have larger disks than local spirals at a given F max , 
whereas values A log < correspond to disks which are smaller than in the 
local universe. As indicated by the fit to the data (solid line), we find a slight 
trend towards smaller disks at higher redshifts, in agreement with theoretical 
predictions within the hierarchical scenario (Mo et al. 1998, dotted line). 



simulated UBgRIJK broad-band colors. Probing a large volume in parame- 
ter space, these synthetic colors were fitted to the observed broad-band colors, 
thus deriving the four model parameters for each of the z > 0.5 FDF galax- 
ies. The best-fitting models indicated that high-mass galaxies on the average 
have higher star formation efficiencies, with a "break" at V max « 140km/s 
which — for reasonable MJL ratios — interestingly is good agreement with 
the results found by Kauffmann et al. (2003) for local galaxies from the Sloan 
Digital Sky Survey. 

Moreover, the best-fitting formation redshifts of the Ferreras et al. models 
were found to be higher for the more massive FDFTF galaxies than for low- 
mass galaxies. The models hence yielded evidence that high-mass spirals 
started to convert their gas into stars at earlier cosmic epochs and on shorter 
timescales than low-mass ones. When evolved to zero redshift, the mean 
model stellar ages turned out to be older in high-mass galaxies than in low- 
mass galaxies. These results hint towards an anti-hierarchical evolution of the 
stellar (baryonic) component, a phenomenon that recently has been referred 
to as "down-sizing" (e.g., Kodama et al. 2004). 
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It is however improbable that this implies a contradiction to the hierar- 
chical growth of the Cold Dark Matter halos. In Fig. El we show the off- 
sets A log rd of the FDFTF galaxies with high quality rotation curves from 
the local velocity-size relation (which we presented in Fig. 01 as a function 
of redshift. Given their maximum rotation velocity, distant galaxies with 
A log rd < have disks that are smaller than in the local universe, while the 
disks of galaxies with A log > are larger. Though the scatter in A log is 
substantial, we observe a slight trend towards smaller disk sizes at higher red- 
shifts. This is in compliance with the results from other observational studies 
(e.g., Giallongo et al. 2000, Ferguson et al. 2004). Moreover, as is depicted in 
Fig. El the fit to the data is in relatively good agreement with the prediction 
of disk growth in the hierarchical scenario (e.g. Mo et al. 1998). 

On the one hand, we find that the observed luminosity evolution of the 
FDF galaxies deviates from the results of simulations which were used to 
predict the TFR of distant spirals. On the other hand, the intermediate- 
redshift disks are observed to be smaller than locally, which is in compliance 
with the CDM hierarchical scenario. The results from single-zone models 
might indicate that the star formation is suppressed in low-mass galaxies due 
to, e.g., SN feedback (cf. Dalcanton et al. 2004). Potentially, the discrepancies 
between observations and simulations arise from the fact that the mechanisms 
suppressing the star formation are not yet implemented realistically enough 
in models of galaxy evolution. 

7 Conclusions 

Using the FORS instruments of the VLT in multi-object spectroscopy mode 
and HST/ACS imaging, we have observed a sample 113 disk galaxies in the 
FORS Deep Field. The galaxies reside at redshifts 0.1 < z < 1.0 and thereby 
probe field galaxy evolution over more than half the age of the universe. 
All spectrophotometric types from Sa to Sdm/Irr are comprised. Spatially 
resolved rotation curves have been extracted and fitted with synthetic velocity 
fields that account for geometric distortions as well as blurring effects arising 
from seeing and optical beam smearing. The intrinsic maximum rotation 
velocities K max were derived for 73 galaxies within the field-of-view of the 
ACS images. Two-dimensional surface brightness profile fits were performed 
to measure the structural parameters like disk inclinations, position angles 
etc. 

The massive distant galaxies fall onto the local Tully-Fisher Relation, 
while the low-mass distant galaxies are brighter than locally by up to >2 m in 
rest-frame B. This trend might be combined with an evolution in metallic- 
ity. We find no evidence for a bias or systematic errors that could induce the 
observed shallow slope of the Tully-Fisher Relation at intermediate redshifts. 
Discrepancies between several previous studies could be explained as a com- 
bination of selection effects and small number statistics on the basis of such 
a mass-dependent luminosity evolution. On the other hand, this evolution 
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would be at variance with the predictions from numerical simulations. For a 
given Vmax, the disks of the distant galaxies are slightly smaller than those 
of their local counterparts, as expected for a hierarchical structure growth. 
Our results therefore are discrepant with theoretical predictions only in terms 
of the stellar populations properties. A possible explanation could be the 
suppression of star formation in low-mass disks which is not yet properly im- 
plemented in models of galaxy evolution. 
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